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代谢 组 与 转录 组 联合 解析 赤 皮 青 内 叶片 黄 化 变异 机 制 
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Abstract: In order to reveal the albino mechanism of gold-coloured mutant leaves of 
Cyclobalanopsis gilva, a naturally-occurring leaf-color mutant was used as experimental materials, 
and the metabolome and transcriptome of mutant leaves and normal green leaves were analyzed 
by UHPLC-Q Exactive HF-X and high-throughput RNA sequencing, respectively. The results 
were as follows: ( 1 ) 257 and 357 significantly changed metabolites ( SCMs ) were respectively 
identified under the positive ion mode and the negative ion mode. Compared with green leaves, 
the content of some flavonoids such as quercetin, leucoyanidin, myricetin and their derivatives 
(pyranodelphinin A, isorhamnetin 3-glucuronide, etc. ) increased significantly in mutant leaves, 
but the content of chlorophylls and carotenoids decreased significantly. ( 2 ) A total of 4 146 
differentially expressed genes (DEGs) were detected, of which 1 711 were up-regulated and 2 435 
were down-regulated. ( 3) KEGG enrichment analysis showed that SCMs and DEGs were mainly 
enriched in photosynthesis, porphyrin and chlorophyll metabolism and flavonoid biosynthesis. 
The research results indicate that inhibition of chlorophyll synthesis, chloroplast developmental 
abnormalities and promotion of flavonoid synthesis were the main factors driving the golden leaf 
coloration in the mutant C. gilva. In addition, the genes of the MYB and bHLH families were 
significantly up-regulated in mutant leaves, confirming these two types of transcription factors 
were involved in regulating flavonoid biosynthesis. This study provided new molecular insights 
into the phenomenon of leaf etiolation, and also provided reference for exploring leaf color-related 
functional genes and breeding of landscape plant. 
Key words: Cyclobalanopsis gilva, leaf-color mutant, albino, metabolome, transcriptome 
叶 色 突 变 是 指 植物 在 生长 发 育 过 程 中 发 生 的 叶 色 变化 现象 〈 刘 新 亮 等 ， 2017) 。 叶 色 突 

变易 于 鉴别 ，Gustafsson 最 早 将 其 分 为 条 纹 、 斑 点 、 浅 绿 、 黄 化 和 白化 $ 类 〈Gustafsson， 
1942) ， 此 后 Manjaya 和 Nandanwar (2007) 又 进一步 细 分 出 黄 绿 、 紫 时、 类 病 斑 等 类 型 。 
黄 化 是 其 中 一 种 重要 的 突变 类 型 ， 黄 化 突变 体 常 被 用 于 植物 光合 生理 《〈 张 晨 劝 等 ，2019) 、 
叶绿体 超 微 结构 〈 李 淑 培 ，2023) 、 叶 绿 素 生 物 合成 Zhu et al.，2014) 等 方面 研究 。 此 外 ， 
黄 化 突变 在 育种 领域 也 有 重要 应 用 , 许多 “ 金 叶 ”类 园林 植物 品种 都 来 自 于 黄 化 突变 , 如 “ 金 
HE Sophora japonica Aurea ) ; “中 华 金 叶 榆 ” 《Ulmus pumila *inye ) ; “ 金 叶 
XS (Acer palmatum *Aurea^) 等 。 目 前 ， 有 关 植 物 叶 色 黄 化 的 研究 集中 于 少数 模式 植 
物 和 一 些 重要 的 农作物 ， 比 如 拟 南 芥 (Arabidopsis thaliana) ; 水 稳 (Oryza sativa) ; 黄瓜 
等 (Li et al.，2012; Maekawa etal., 2015; 能 兴 伟 等 ，2023) 。 在 园林 植物 中 ， 仅 在 黄山 
Ab] (Koelreuteria bipinnata var. integrifoliola) (Lyu etal,2017) ; {RA (Ginkgo biloba) 

(Li etal., 2019) ; 杂交 构 树 (Broussonetia kazinokiX B. papyrifera) (Wang et al., 2022) 等 
少数 植物 中 有 过 报道 。 
近年 来 , 高 通 量 测序 技术 被 广泛 用 于 生物 组 学 研究 , 能 够 从 分 子 层面 揭示 生物 现象 和 生 
物 过 程 的 发 生机 制 ( 熊 兴 伟 等 ， 2023) 。 并 且 ， 通 过 多 组 学 的 联合 应 用 ， 可 以 更 加 直观 地 反 
映 机 体 的 变化 水 平 ,为 生物 体 静 态 和 动态 的 变化 提供 更 为 深入 和 广阔 的 研究 视野 ( 陆 小 雨 等 ， 
20200 . Lyu 等 (2017) 通过 高 通 量 测序 技术 ， 从 黄山 栾 树 品种 “ 金 焰 彩 栾 ” CK. bipinnata 
var. Integrifoliola ‘Jinyan’) 中 检测 出 9 个 叶绿素 代谢 和 14 个 类 胡萝卜 素 生物 合成 的 关键 基 
Kle Wang 等 〈2022) 采用 多 组 学 方法 对 杂交 构 树 的 黄 叶 突变 体 进行 研究 ， 揭 示 了 金黄 叶 色 
表 型 的 形成 与 光合 色素 的 含量 与 比例 变化 以 及 叶绿体 结构 和 功能 缺陷 有 关 。Yamashita 等 
(2021) 对 茶树 黄 化 品种 (Camellia sinensis Koganemidori ) 的 叶绿素 代谢 和 和 氨基酸 积 累 机 
制 进行 了 探究 ,表明 供 试 品种 的 黄 化 表 型 是 由 于 叶绿体 发 育 和 叶绿素 合成 相关 基因 缺乏 所 致 ， 
高 氨基 酸 则 源 于 代谢 相关 基因 的 上 调 表 达 所 引起 的 广泛 和 蛋白质 降解 。Luo 等 (20220 采用 转 


mh 
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录 组 和 代谢 组 联 用 技术 , 揭示 了 甘蔗 叶 色 黄 化 现象 与 叶绿素 合成 途径 减少 、 光 合 基因 表达 下 
降 、 金 属 离子 调控 功能 障碍 以 及 次 生 代谢 物质 改变 等 因素 有 关 。 总 而 言 之 ,植物 时 色 黄 化 与 
叶绿素 代谢 相关 基因 的 突变 或 表达 受阻 有 重要 关系 。 此 外 ,也 会 受到 花 青 素 、 类 胡萝卜 素 等 
其 他 色素 代谢 变化 的 综合 作用 〔〈 林 馨 颖 等 ，2022) 。 

赤 皮 青冈 (Cyclobalanopsis gilva) 为 过 斗 科 (Fagaceae) 青冈 属 CCyclobalanopsis) 党 
绿 乔 木 , 主要 分 布 于 我 国 南方 海拔 300~1 500 m 的 山林 地 带 (中 国 植物 志 编 辑 委员 会 , 1998). 
其 木材 质地 坚硬 ， 纹 理 细 腻 ， 是 优良 的 硬木 用 材 树种 ， 也 可 用 于 低 山 丘陵 混交 造林 、 林 下 补 
植 和 园林 绿化 (欧阳 泽 怡 等 ，2021; 秦 之 旷 等 ，2023) 。 本 研究 团队 于 2016 年 在 宁波 海曙 
溪 下 育苗 基地 发 现 一 株 赤 皮 青 内 播种 变异 株 , 其 成 熟 叶片 为 黄色 , 新 叶 金 黄色 , 枝 干 明黄 色 ， 
并 于 2018 年 底 通 过 嫁接 繁殖 多 株 材 料 ， 其 特异 性 状 已 连续 4 年 表现 稳定 。 本 研究 以 赤 皮 青 
闷 黄 化 突变 植株 和 正常 植株 的 叶片 为 研究 材料 ,采用 代谢 组 和 转录 学 联 用 技术 同时 结合 生理 
生化 测定 ， 拟 探讨 以 下 问题 : (1) 赤 皮 青冈 变异 植株 金黄 叶 色 表 型 形成 的 生理 机 制 ， (2) 
变异 植株 叶 色 黄 化 的 转录 调控 机 制 。 本 研究 可 为 金 叶 系 园林 植物 的 选 育 和 遗传 改良 提供 参考 
依据 。 


1 材料 与 方法 


1.1 试验 材料 

以 赤 皮 青冈 黄 化 突变 体 (Yellow leaf, YL) 和 正常 植株 (No yellow leaf, NYL) 为 研究 
材料 ， 两 种 材料 均 种 植 于 宁波 市 农业 技术 推广 总 站 林 特 基地 C121242/24"E, 29948'46"ND 。 
2021 年 8 月 采集 两 种 植株 叶片 样品 (图 1) ， 先 置 于 液 氮 速冻 ， 后 置 于 -80C 冰 箱 中 保存 备 
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图 1 赤 皮 青冈 黄 化 植株 YL(A) 和 叶 色 对 比 图 (B) 
Fig. 1 Gold-coloured mutant YL ( A ) and leaf color contrast ( B ) 

1.2 代谢 组 分 析 
1.2.1 代谢 物 提取 

称 取 50 mg 样品 至 离心 管 ， 加 入 400 uL 的 甲醇 : 乙 且 (VV=1:1， 含 内 标 ) 提取 液 。 
用 冷冻 组 织 研磨 仪 (Wonbio-96c， 上 海 万 柏 ) 进行 研磨 ， 提 取 液 静 置 后 在 4C 条 件 下 13 000 
g 离心 15 min, 之 后 取 上 清 液 进行 上 机 分 析 。YL M NYL 组 各 设置 6 个 生物 学 重复 。 每 个 
样品 取 20 uL 上 清 液 ， 混 合 后 用 于 质 控 分 析 。 
1.2.2 代谢 物 检 测 

利用 Thermo Fisher Scientific 超 高 效 液 相 色谱 (CUHPLC ) 串 联 侍 里 叶 变 换 质 谱 (Q Exactive 
HF-XO 进行 LC-MS 分 析 。 色 谱 柱 : ACQUITY UPLC HSS T3 (100 mm x 2.1 mm, 1.8 um; 
Waters, USA) 。 液 相 条 件 : 流动 相 A 为 95% 水 -59eZ.H CE 0.195 P EO ,流动 相 B 为 


47.5% 乙 且 + 47.5% 异 丙 醇 + 5% 水 〈 含 0.1% 甲 酸 ) 。3) 洗 脱 梯度 : 正 离子 模式 : 0 min, 
A/B (V/V) 为 100:0; 3.0 min, 为 80:20; 4.5 min, 为 65:35; 5.0 min, 为 0:100; 6.3 min, 
X 0:100; 6.4 min, 为 100 :0; 8.0min， 为 100 :0。 负 离子 模式 : 0min, A/B (V/V) 为 
100:0; L5 min, 为 95:5; 2.0 min, X 90:10; 4.5 min, X 70:30; 5.0 min, 730:100; 
6.3 min, 73 100:0; 6.4 min, 73 100:0; 8.0 min, X 100:0. 4) 柱 温 40'C; 进 样 量 3 uL; 
流速 0.4 mL-min-!。 质谱 条 件 : 电 喷雾 离子 源 CESD ， 喷 雾 电 压 3 500 V OESR) 和 -3 500 
V〔 负 模式 ) ; 加 热 温度 425C; 毛细 管 温度 325C; WAE 50 Arb; 辅助 气流 速 13 Arb; 
扫描 范围 70 ~ 1 050 mz; 一 级 分 辨 率 60 000; 二 级 分 辩 率 7 500; 碰撞 能 分 别 为 20、40、 
60 eV; 采用 质谱 连续 扫描 采集 数据 。 

1.2.3 数据 分 析 

将 原始 数据 导入 Progenesis v2.2 软件 (Waters，USA ) ， 通 过 与 HMDB、METLIN 等 代 

谢 数据 库 比 进行 代谢 物 鉴定 。 利 用 R 软件 包 进 行 代 谢 组 数据 分 析 和 热 图 制作 ， 对 两 组 样品 
进行 主 成 分 分 析 Cprincipal component analysis; PCA ) 和 正 交 偏 最 小 二 乘法 分 析 (OPLS-DA)， 
得 到 模型 的 变量 权重 值 VIP Cvariable important in projection? , LA VIP > 1, P<0.05 并 且 FC 

(Fold Change) > 1.10 或 FC < 0.91 为 标准 盘 选 组 间 的 显著 差异 代谢 物 (significantly changed 
metabolites, SCMs) ， 通 过 与 KEGG 数据 库 进行 代谢 途径 富 集 。 
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1.3 转录 组 测序 与 数据 处 理 
1.3.1 转录 组 测序 及 文库 构建 

采集 YL 和 NYL 的 新 鲜 叶 片 ， 参 考 陆 小 雨 等 (20200 方法 进行 RNA 提取 和 质 控 处 理 。 
采用 TruseqIMRNA sample prep Kit (Illumina ) 的 方法 构建 文库 。 首 先 ， 用 带 Oligo (dT) 的 
磁 珠 分 离 mRNA。 接 着 ,将 mRNA 随机 打 断 ,在 SMART 逆转 录 酶 (SMARTScribeTM Reverse 
Transcriptase) 作用 下 合成 cDNA 第 1 条 链 ， 再 利用 RNaseH 将 RNA 链 降 解 ， 合 成 cDNA 
第 2 条 链 。 之 后 ， 通 过 End-Repair Mix 将 双 链 cDNA 补 齐 ， 其 3” 末端 加 A 处 理 后 连接 测 
序 接头 。 然 后 ,将 cDNA 进行 PCR. 富 集 ,再 利用 AMPure XP beads 纯化 产物 。 在 经 过 TBS-380 

(Picogreen) 定量 后 ， 按 数据 比例 混合 上 机 。 测 序 时 先 通过 cBot 平台 完成 桥 式 PCR 扩 增 ， 

生成 艇 (clusters )， 再 利用 Ilumina Noveseq 6000 平台 完成 转录 组 测序 ， 构 建 Iumina PE X 
库 〈 读 长 2bpx150bp) o 
1.3.2 数据 处 理 

将 测序 得 到 的 图 像 信号 经 CASAVA 碱 基 识 别 转换 为 原始 数据 Craw data) ， 经 过 质 控 
获得 到 高 质量 的 干净 数据 (clean data) 。 采 用 无 参考 基因 组 的 转录 组 分 析 , 通过 Trinity v2.8.5 
软件 (Grabherr et al., 2011) 将 Clean reads E22 EE IE. Contig) 和 单一 序列 〈Singleton ) 。 
此 后 ， 利 用 TransRate v1.0.3 (Smith-unna et al., 2016) 和 CD-hit v4.5.7 (Li & Godzik, 2006) 
对 初始 组 装 序 列 进 行 优化 过 滤 ， 并 利用 BUSCO v3.02 {iF (Simão et al., 2015) 再 次 评估 ， 
得 到 后 续 分 析 的 最 终 转 录 本 。 
将 转录 本 与 6 大 数据 库 (NR, Swiss-Prot, Pfam, eggNOG, GO M KEGG) 比 对 以 获 
取 注 释 信 息 。 通过 RSEM v1.3.1 软件 (Li et al., 2011) 计算 基因 的 TPM (transcripts per million 
reads) 表达 量 ， 利 用 DESeq2 v1.24.0 软件 (Michael et al., 2014) 筛选 2 个 组 的 差异 表达 基 
因 (differentially expressed genes, DEGs) ， 将 阀 值 差异 倍数 FC >2 且 修 正己 值 CP adjust 
< 0.05 作为 筛选 标准 。 以 P adjust < 0.05 为 标准 , 采用 Goatools v0.6.5 软件 (Klopfenstein et al., 
2018) 进行 GO 富 集 分 析 和 KEGG pathway 富 集 分 析 ， 筛 选 参与 叶 色 变 化 的 关键 基因 。 


1.3.3 代谢 物 与 基因 相关 性 分 析 
采用 Prism 8.0 (GraphPad, USA) 软件 进行 SCMs 与 DEGs 的 皮尔 森 (pearson) 相关 


性 分 析 。 
关系 数据 。 
1.3.4 qRT-PCR 验证 


以 相关 性 系数 |r| > 0.8 HP «0.05 为 阔 值 ， 得 到 SCMs 与 DEGs 的 分 子 间 相互 作用 


选用 8 个 DEGs 进行 RT-qPCR 来 检验 转录 组 数据 的 可 靠 性 ， 利 用 Primer Premier 5.0 


软件 设计 PCR 引物 ( 表 1) o EA C4Cs (Accession number: ID67285000 为 内 参 基 因 


(Marum 


et al., 2012) 。 设 置 3 次 生物 学 习 


的 


复 。 将 荧光 定量 表达 水 平 进行 Log: 转换 ， 与 对 应 基 医 


转录 组 表达 丰 度 进行 比较 。 
表 1 qPCR 引物 序列 
Table 1 Primer sequences for qPCR 
pj 。 正 疝 引物 序列 (3) 反问 引物 序列 C5 23» 
Gene ID Forward primer sequence (5 一 3 ) Reverse primer sequence (5 一 3 ) 
CACs TCTGGGAGAAGAGTGGCTACA GAGCCACCATTCAAATCCT 
FLS CCTCCTTCCACCATCAACTA ATGACCTTCAACCCCTAACC 
CHS CAGCAACAAAGGCAATCAA GACAAGAACACGAGCACCC 
UROD CTTCTTTCTTTGCTCCTTT GGTCTTTCAACATCTTCCC 
psaL GTCCAATCTGAAAAGCCAA GTCCAACAAGGAGAAACCC 
psaV CTCTTCCTCCACCTTGTTC GTGCTTGACTTGGTTTTTG 
ELIP CAACACCTCAGCCACCACA GCAGCAACAAAGCCAATCA 
MYB7 TCACCTGCCAAAGCCCGTT CCGACCCCATTTTCACCGT 
Cab GGAACTAAAGAATGGGAGG CAGCAATATGGTCGAAGAG 
1.4 叶绿素 含量 测定 


取 YL 和 NYL 的 新 鲜 叶 片 ， 经 过 
来 测定 两 组 叶片 中 的 叶绿素 含量 。 


1.5 光合 参数 测定 


采用 CI-340 便携 式 光 合作 


(omg. f 


8. 
2 结果 与 分 析 


2.1 代谢 组 测序 结果 与 分 析 
2.1.1 代谢 组 的 多 元 统计 分 析 


对 YL 和 NYL 两 组 样本 


M 


物 


ETE 


、 负 


离子 模式 


磨 和 过 滤 后 ， 参 照 张 丽 霞 等 (2021) 的 方法 
设置 3 次 生物 学 重复 。 


测定 仪 (CID, USA) 对 YL 突变 体 和 NYL 植株 叶片 的 净 
光合 作用 、 气 孔 导 度 、 胞 间 CO 浓度、 蒸腾 速率 等 光合 参数 进行 测定 。 设 置 3 次 生物 学 重 


下 的 UHPLC-Q Exactive HF-X 数据 进行 代谢 
成 的 多 元 统计 分 析 。PCA 分 析 结 果 显 示 ， 正 离子 模式 下 第 1 ERI (PC1) 和 第 2 ER 


分 (PC2) 分 别 占 总 变量 的 40.60% 和 19.10%， 负 离子 模式 下 第 1 主 成 分 (PC1) 和 第 2 E 
2: A，B) 中 。OPLS-DA 分 析 结 果 显 示 


成 分 (PC2 ) 分 别 占 总 变量 的 37.30% 和 24.10% (B 
2: C, DO ， 两 种 离子 模式 下 差异 代谢 物 R2 


(图 
0， 说 明 样 本 数据 描述 恨 好。 


在 正 、 负 离子 模式 下 


F) 。 


鉴定 出 614 个 SCMs, H: I 
个 ， 下 调 的 有 109 个 ; 负离子 模式 下 357 个 ， 上 调 的 有 147 个 ， 下 调 的 有 210^ CER 


直 均 高 于 Q2 值 ， 且 Q2 5 Y f 


和 截 距 小 于 
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QC. Quality control; VIP. Variable importance in projection. 


图 2 黄 化 植株 YL 和 正常 植株 NYL 的 PCA 得 分 图 (A、B)、OPLS-DA 得 分 图 (C、D) 和 差异 


代谢 物 火山 图 (E，F) 
Fig.2 PCA score plots ( A, B ) OPLS-DA scores plots ( C, D )and volcano plots ( E, F ) of 
differential metabolites between YL and NYL 
2.1.2 主要 的 SCMs 分 析 
从 614 个 SCMs 中 筛选 出 相对 含量 前 30 差异 代谢 物 〈 表 2) ， 其 中 : 黄酮 类 化 合 物 有 9 
种 ， 均 在 YL HARA EH, ARE A 上 调 最 显著 ， 上 调 2.28 fü. MER 3-0- 〈6"- 乙 
酰 葡 萄 糖苷 ) 也 上 调 了 1.55 倍 ; 5 种 核 音 及 其 衍生 物 在 YL 中 有 差异 积累 ， 其 中 2 种 喀 喧 核 
Tux EU, 3 种 嗓 叭 核 痛 显著 下 调 ; 9 种 脂 类 代谢 物 中 ， 有 5 种 在 YL 中 受到 上 调 ， 主 要 
为 异 戊 烯 醇 脂 类 ; 2 种 有 机 酸 中 ，1 种 显著 上 调 ， 另 1 种 显著 下 调 。 


表 2 YL 和 NYL 组 中 相对 含量 前 30 的 SCMs 
Table 2 The top 30 SCMs between YL and NYL 


nm m NL ZR DLOSL um 
吡 喃 酮 啡 肽 A Pyranodelphinin A TKE Flavonoids C30H33016 1.45E-05 3.541 1 2.280 8 上 调 Up 
HER 3-O- 咖 啡 酰 和 葡萄 糖苷 Quercetin 3-O-caffeyl-glucoside iKAW Flavonoids C30H26015 4.24E-06 3.304 1 1.522 3 上 调 Up 
麦迪 考 赛 德 G Medicoside G 异 戊 烯 醇 脂 Prenol lipids CaligDis 00398 31865 17847 ”上调 Up 
SIME Isopentenyl adenosine In EF Purine nucleosides CisHaiNsO4 0.003 6 3.1050 0.513 ”下 调 Down 
MRE Oleragenoside 异 戊 烯 醇 脂 Prenol lipids CeHa4Oi6 0.037 2 2.858 2 1.535 8 上 调 Up 
WERE Deoxyuridine REKET Pyrimidine nucleosides CoHi2N205 0.000 9 2.785 7 1.407 7 上 调 Up 
乳酸 菌 素 Lactucin 内 酯 类 Lactones CisHi6Os 0.002 1 2.748 4 0.5605 ”下 调 Down 
20- 氧 代 Kis E4 20-Oxo-leukotriene E4 站 肪 醇 Fatty Acyls CzH35NO6S 0.012 8 2.7385 1.6553 上 调 Up 
乳糖 -N- 二 糖 工 Lacto-N-biose I 有 机 氧化 合 物 Organooxygen compounds C14H25NO11 0.0145 2.695 2 1.505 4 上 调 Up 
原 花 青 素 Proanthocyanidins KAW Flavonoids C31H28012 0.008 0 2.683 3 1.439 3 上 调 Up 
13E- 四 聚 -16- 羧 基 -LTE4 13E-Tetranor-16-carboxy-LTE4 羧 酸 及 衍生 物 Carboxylic acids and derivatives CioH27NO7S 0.000 6 2.639 3 1.404 2 上 调 Up 
西 红 花 苷 IV Dicrocin 异 戊 燃 醇 脂 Prenol lipids CoHaOi4 0.009 9 2.6181 0.6970 MÑ Down 
HEZA 3-0-(6"- 乙 酰 葡萄 糖苷 ) Quercetin 3-O-(6"-acetyl-glucoside) ”黄酮 类 化 合 物 Flavonoids C5H2015 0.006 1 2.5950 1.547 9 上 调 Up 
T- : 小 其 5 
ive uem Ae 羧 酸 及 衍生 物 Carboxylic acids and derivatives CHHI7NOs 6.90E-09 2.582 7 0.824 1 下 调 Down 


2-Deoxy-2,3-dehydro-n-acetyl-neuraminic acid 


苯 甲 酰 芍 药 苷 Benzoylpaeoniflorin 异 成 烯 醇 脂 Prenol lipids C30H32012 0.002 8 2.5750 1.359 6 上 调 Up 
燕 草 素 -3-0- 葡 萄 糖 基 -5-O0- 咖 啡 栈 葡 萄 糖苷 3 : 
黄酮 类 化 合 物 Flavonoids C36H36020 9.24E-05 2.564 7 1.259 5 上 调 Up 
Delphinidin 3-O-glucosyl-5-O-caffeoylglucoside 
AREJA Capillarisin 香 豆 素 及 其 衍生 物 Coumarins and derivatives Ci6H1207 4.85E-06 2.560 3 1.307 9 上 调 Up 
阿 萨 米 卡 因 A Assamicain A 黄酮 类 化 合 物 Flavonoids Ca4H36022 0.000 2 2.5529 1.262 0 上 调 Up 
溶血 磷脂 酰 肌 醇 Lysophosphatidylinositol 甘油 磷脂 Glycerophospholipids CHH2IOIDP 0.0007 2.544 6 1.290 5 上 调 Up 
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2.1.3 SCMs 的 KEGG 通路 分 析 
将 YL 和 NYL 的 SCMs 进行 KEGG 富 集 分 析 ， 共 富 集 到 66 条 通路 ， 主 要 包含 甘油 磷 
脂 代谢 、 辅 助 因子 生物 合成 、 类 黄酮 生物 合成 、ABC 转运 蛋白 、 异 黄酮 生物 合成 、 氨 基 糖 
和 核 苷 酸 糖 代 谢 等 13 条 代谢 通路 〈 图 3) 。 其 中 ， 黄 酮 类 生物 合成 富 集 的 SCMs 最 多 ， 包 
括 类 黄酮 生物 合成 富 集 12 个 , 异 黄酮 生物 合成 富 集 7 个 , 黄酮 和 黄酮 醇 生物 合成 富 集 4 个， 
暗示 黄酮 类 物质 可 能 与 YL 叶 色 黄 化 有 关 。 此 外 ， 甘 油 磷脂 代谢 、 辅 酶 合成 代谢 、ABC 转 
运 和 蛋白 、 氮 基 糖 和 核 苷 酸 糖 代 谢 以 及 莱 丙 烷 生 物 合成 等 分 别 富 集 了 17 个、12 个 、7 个 、56 
个 和 4 个 代谢 物 。 


类 黄酮 生物 合成 Flavonoid biosynthesis] Ge P 值 
甘油 磷脂 代谢 Glycerophospholipid metabolism1 e Pvalue 
异 黄酮 生物 合成 Isoflavonoid biosynthesis 4 e 0.040 
Mi uie e EL ; i 0.030 
植物 激素 信号 转 导 Plant hormone signal transduction -4 © 0.020 
核 黄 素 代 谢 Riboflavin metabolism 4 e f: 
ABC 转运 蛋白 ABC transporters] — 6 0 
黄酮 和 黄酮 醇 生 物 合成 Flavone and flavonol biosynthesis 4 o 
丙 氨 酸 、 天 冬 氨 酸 和 谷 氨 酸 代谢 Alanine, aspartate and glutamate metabolism 4 e 
AMAZE ARI Amino sugar and nucleotide sugar metabolism] @ 数量 Number 
辅助 因子 生物 合成 Biosynthesis of cofactors 4 9 e 5.25 
烟 酸 和 烟 酰 胺 代谢 Nicotinate and nicotinamide metabolism ] e e 7.50 
玉米 素 生物 合成 Zeatin biosynthesis 4 o e ^75 
葵 两 烷 生 物 合成 Phenylpropanoid biosynthesis 4 [z e 0 


0 0.05 0.10 0.15 0.20 0.25 0.30 
富 集 因子 Rich factor 
图 3 黄 化 植株 YL 和 正常 植株 NYL 差异 代谢 物 的 富 集 图 
Fig.3 Enrichment of different metabolites between YL and NYL 

2.2 转录 组 测序 结果 与 分 析 
2.2.1 测序 质量 分 析 

对 NYL 和 YL 材料 的 RNA 进行 转录 组 测序 ， 质 控 后 共 获 得 4426 Gb 干净 数据 (Clean 
data) ( 表 3) 。 各 样品 的 干净 序列 (Clean read) 条 数 在 43 864 122 和 56 594 944 之 间 ，Q20 
直 介 于 97.61% 和 97.79%, Q30 值 介 于 93.03% 和 93.49%, GC 碱 基 占 比 在 44.49% 和 44.99% 
之 间 ， 满 足 于 后 续 分 析 要 求 。 


-— 


表 3 转录 组 测序 质量 分 析 


Table 3 Quality analysis of transcriptome sequencing of Cyclobalanopsis gilva 


干净 数据 错误 率 
样品 名 称 干净 序列 Q20 比例 Q30 比例 GC 含量 
Clean data Error rate 
Sample ID Clean reads Q20 ratio (%) Q30ratio (%) GC content (%) 
(bp) (%) 
YL-1 6 712 458 677 45 836 326 0.03 97.61 93.03 44.50 
YL-2 7 689 229 888 52 647 488 0.03 97.69 93.31 44.99 
YL-3 7 372 128 382 50 375 710 0.03 97.64 93.07 44.49 
NYL-1 6 381 564 905 43 864 122 0.03 97.70 93.25 44.73 
NYL-2 8 227 812 354 56 594 944 0.03 97.74 93.31 44.76 
NYL-3 7 871 942 680 53 811 634 0.03 97.79 93.49 44.78 
2.2.2 差异 表达 基因 统计 


从 YL 和 NYL 的 样品 中 分 别 鉴定 了 46 391 和 48 018 个 表达 基因 ,其 中 16 699 和 15 072 
个 基因 分 别 在 YL M NYL 中 特异 表达 (图 4: A) 。 以 NYL IEH, Æ YL 中 共 检 测 到 4 
146 个 DEGs， 其 中 1711 个 上 调 ，2 435 个 下 调 (图 4: BD 。 


2 435 
B 上 调 Up 


B Mi Down 


Number of DEGs 
T 
© 
© 
L 


EDGs. 差异 表达 基因 ; Up. 上 调 ; Down. 下 调 。 
EDGs. Differentially expressed genes; Up. Up-regulated; Down. Down-regulated. 
图 4 YLA NYL 组 间 差 异 表 达 基 因 的 韦 恩 图 (A) 和 差异 表达 基因 上 /下 调 统计 图 (B) 
Fig. 4 Venn diagram ( A ) and differentially expressed gene up/down regulation ( B ) between 
NYL and YL 


2.2.3 DEGs 的 GO 功能 分 析 

GO 分 析 表 明 ，YL fü NYL 的 3 612 个 DEGs (P<0.05) 显著 富 集 到 12 个 生物 学 过 程 、 
14 个 细胞 组 分 和 2 个 分 子 功能 条 目 〈 图 5) 。 在 生物 学 过 程 中 ， 富 集 基因 最 多 的 是 “前 体 代 
谢 产物 和 能 量 产生 ”， 其 次 为 “光合 作用 ”、“ 色 素 生物 合成 过 程 ”、“ 色 素 代谢 过 程 ”、 
“中 啉 化 合 物 生物 合成 过 程 ” 等 ,涉及 过 程 多 与 光合 作用 或 叶绿素 合成 有 关 ， 以 下 调 基 因为 
主 ， 也 有 少数 与 磷酸 脱氧 酶 、 质 体 丙 酮 酸 激酶 、 葡 萄 糖 基 转 移 酶 及 二 氧 黄酮 醇 还 原 酶 相关 的 
代谢 途径 发 生 了 不 同 程度 上 调 。 在 细胞 组 分 中 , 富 集 最 多 的 GO 条 目 是 “ 质 体 ”和 “叶绿体 ”， 
其 次 为 “类 宫 体 膜 ”“ 光 合 膜 ”“ 叶 绿 体 类 宫 体 膜 ” 等 ， 也 以 下 调 基 因 的 数目 较 多 。 在 分 子 
功能 中 ， 存 在 显著 差异 的 代谢 途径 最 少 ， 仅 有 “四 吡咯 结合 ”和 “叶绿素 结合 ”， 都 以 下 调 
基因 居多 ， 其 中 “叶绿素 结合 ”途径 当中 的 下 调 基因 占 总 数 的 90% 以 上 。GO 功能 富 集 分 析 
KRH, YL 的 叶 色 黄 化 突变 与 光合 作用 、 叶 绿 素 代 谢 等 过 程 有 密切 关系 。 
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l. 前 体 代谢 产 物 和 能 量 产生 ; 2. 光合 作用 ; 3. 色素 生物 合成 过 程 ; 4. 色素 代谢 过 程 ; $. 中 啉 化 合 物 生 物 


合成 过 程 ; 6. 四 吡咯 生物 合成 过 程 ; 7. 光合 作用 , 光 捕 获 ，8. 光合 作用 , 光 系 统 I 光 捕 获 ; 9. 叶绿素 生物 合 
成 过 程 ; 10. 叶绿素 代谢 过 程 ; 11. 细胞 对 光 刺 激 的 反应 ; 12. 光合 系统 I 中 光合 电子 传输 ; 13. 质 体 ; 14. 叶 
绿 体 ，15. 类 圳 体 膜 ，16. 光合 膜 ，17. 叶绿体 类 圳 体 膜 ，18. 质粒 类 宫 体 膜 ，19. 光 系 统 ，20. 光 系 统 I; 
21. 光 系 统 I: 22. 光 系 统 工 放 氧 复合 体 ; 23. 质 体 小 球 ; 24. 光 系统 I 反应 中 心 ; 25. 质 体 类 圳 体腔 ; 26. 叶 
RREK: 27. 四 吡咯 结合 ，28. 叶绿素 结合 。 


1. Generation of precursor metabolites and energ; 2. Photosynthesis; 3. Pigment biosynthetic process; 4. Pigment 
metabolic process; 5. Porphyrin-containing compound biosynthetic process; 6. Tetrapyrrole biosynthetic process; 
7. Photosynthesis, light harvesting ; 8. Photosynthesis, light harvesting in photosystem I; 9. Chlorophyll 
biosynthetic process; 10. Chlorophyll metabolic process; 11. Cellular response to light stimulus; 12. Photosynthetic 
electron transport in photosystem I; 13. Plastid; 14. Chloroplast; 15. Thylakoid membrane; 16. Photosynthetic 
membrane; 17. Chloroplast thylakoid membrane; 18. Plastid thylakoid membrane ; 19. Photosystem ; 20. 
Photosystem I; 21. Photosystem II; 22. Photosystem II oxygen evolving complex; 23. Plastoglobule; 24. 
Photosystem I reaction center; 25. Plastid thylakoid lumen; 26. Chloroplast thylakoid lumen; 27. Tetrapyrrole 
binding; 28. Chlorophyll binding. 
图 5 YL 组 和 YL 组 差异 表达 基因 GO 富 集 分 析 
Fig. 5 GO enrichment analysis on differentially expressed genes in YL and NYL 

2.2.4 差异 代谢 途径 的 KEGG 富 集 分 析 

KEGG 富 集 分 析 结 果 显 示 ， 共 有 744 条 DEGs 被 注释 到 KEGG 通路 ， 其 中 192 条 基因 
显著 富 集 到 植物 -病原 体 相 互 作 用 、 光 合作 用 、 乙 醛 酸 及 二 羧 酸 代 谢 等 12 个 过 程 ( 表 4) 。 
与 光合 作用 相关 的 途径 有 5 条 ， 其 中 : 光合 作用 -天 线 蛋 白质 富 集 因子 最 高 ， 为 0.666 7， 该 
通路 注释 到 10 个 基因 ， 全 部 为 捕 光 叶 绿 素 蛋白 复合 体 〈light-harvesting chlorophyll protein 
complex, LHC) 编码 基因 ; 富 集 到 光合 作用 的 基因 有 23 个 ， 主 要 涉及 光 系 统 的 反应 中 心 复 
AH (reaction-center complex, RCC) ; 富 集 到 光合 生物 碳 固定 途径 的 有 16 个 ， 包 括 磷酸 
丙 糖 异 构 酶 〈triosephosphate isomerase, TIM) . 3-WER?H W R Mi NE (glyceraldehyde 
3-phosphate dehydrogenase, GAPA) 等 关键 酶 的 编码 基因 ; 富 集 到 中 啉 和 叶绿素 代谢 的 有 
12 个 ; 富 集 到 乙 醛 酸 及 二 羧 酸 代谢 的 则 有 22 个 基因 。 此 外 ， 一 些 DEGs 也 参与 了 类 黄酮 生 
物 合 成 、 花 青 素 生物 合成 、 氮 基 酸 代谢 、 葵 丙烷 生物 合成 等 途径 。 


KA YL 和 NYL 间 显 著 富 集 KEGG 通路 (修正 P< 0.05) 
Table 4 Significantly enriched KEGG pathways between YL and NYL ( P adjust < 0.05 ) 


"5 


Zps ra B Jo "8 "8 M > 大 其 通路 注释 的 所 5 5 内 = 4E RE VE TE fd 
编写 通路 名 称 通路 ID 通路 注释 的 差异 表达 基因 AI enes with pathway IT P (B 尸 修正 什 
No. Pathway Pathway ID DEGs with pathway annotation A P Y Rich factor P value P adjust 
annotation 
1 光合 作用 Photosynthesis map00195 23 55 0.418 2 7.01E-15 7.71E-13 
2 光合 作用 一 天 线 蛋 白质 Photosynthesis-antenna proteins map00196 10 15 0.6667 8.61E-10 4.73E-08 
3 植物 -病原 体 相互 作用 Plant-pathogen interaction map04626 62 545 0.113 8 1.79E-07 6.55E-06 
4 hp E SERIE Porphyrin and chlorophyll metabolism map00860 12 60 0.200 0 0.000 130 233 — 0.003 581 409 
5 乙 醛 酸 及 二 羧 酸 代谢 Glyoxylate and dicarboxylate metabolism map00630 22 193 0.1140 0.001 664 845 — 0.030 522 151 
光合 生物 碳 固定 途径 Carbon fixation in photosynthetic 
0 map00710 16 121 0.1322 0.001 515 940 — 0.033 350 685 
organisms 
7 AE tRNA 生物 合成 Aminoacyl-tRNA biosynthesis map00970 17 140 0.121 4 0.002 780 452 0.043 692 822 
8 RWE WX Phenylpropanoid biosynthesis map00940 11 191 0.0576 0.541 863 137 0.961 370 081 
9 花 青 素 生 物 合成 Anthocyanin biosynthesis map00942 2 2 1.000 0 0.003 300317 . 0.045 379 358 
10 ”类 黄酮 生物 合成 Flavonoid biosynthesis map00941 10 66 0.151 5 0.004 205 476 — 0.046 260 234 
HAR, ARMI ARH Glycine, serine and threonine 
11 map00260 15 121 0.124 0 0.003 944 472 0.048 210 213 
metabolism 
12 ”黄酮 和 黄酮 醇 生物 合成 Flavone and flavonol biosynthesis map00944 3 7 0.428 6 0.004 561 683 — 0.045 616 834 


2.2.5 差异 转录 因子 鉴定 与 分 析 

转录 因子 因 与 功能 基因 调控 区 域 结合 而 影响 基因 表达 , 会 对 许多 生物 学 过 程 产生 影响 Cl] 
飞 等 ，2022) 。 对 4 146 个 DEGs 中 的 转录 因子 基因 进行 鉴定 ， 共 筛 得 39 个 转录 因子 相关 
基因 ， 属 于 15 个 家 族 (图 6) 。 其 中 ， 差 异 基因 数目 最 多 的 是 MYB 家 族 (8 个 DEGs) ， 
其 次 为 AP2/ERF Ji bHLH, 分别 有 5 个 和 4 个 ，WRKY、bZIP 和 SBP 则 都 有 3 个 。15 个 
转录 因子 家 族 中 , WRKY、GRAS、AP2/ERF 和 C2H2 以 下 调 基因 居多 ， 其 他 则 以 上 调 为 主 。 
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差异 基因 数目 Number of DEGs 


图 6 差异 基因 中 不 同 转录 因子 统计 分 析 


Fig.6 Statistical analysis of transcription factors in DEGs 


2.2.6 DEGs 的 RT-qPCR 验证 分 析 
选取 8 个 DEGs 进行 qPCR 验证 (图 7) 。 结 果 表 明 ， 痰 光 定 量 表达 与 转录 组 基因 表达 
的 变化 趋势 一 致 ， 证 实 转录 组 数据 可 靠 。 


8 EGG q-PCR — — RNA-seq 8 
S 2 
T vE 
33 g RF 
RG ELE 
RL m Dd 
m-s Ew u 
上 总 2g 
5p pr 
& ed HE 
E "HS 
& "M 2 
总 E. i S 
总 E- 
FLS CHS UROD  psaL psaV ELIP MYB7 Cab 
7 YL 和 NYL 中 8 个 DEGs 的 表达 水 平 验证 
Fig.7 Verification of expression levels of 8 DEGs in YL and NYL 
2.3 YL 的 DEGs 和 SCMs 关联 分 析 


2.3.1 色素 代谢 相关 DEGs 和 SCMs 的 关联 分 析 

植物 叶片 颜色 主要 由 叶绿素 、 类 胡 葛 卜 素 、 花 青 素 等 植物 色素 决定 (林志颖 等 ，2022 )。 
通过 测定 赤 皮 青冈 YL 和 NYL 的 叶绿素 和 类 胡 葛 让 素 含量 , 发 现 YL 中 叶绿素 a(0.35 mg g D. 
叶绿素 b (0.31 mgg!) 、 总 叶绿素 (0.66 mgg!) URAH h% (0.05 mgg!) 的 含量 都 


较 NYL (分 别 为 2.38、1.49、3.87 和 0.44 mgg!) 显著 下 降 。 此 外 ， 叶 绿 素 a/b 的 比值 也 显 
著 降 低 〈 图 8) ， 表 明 YL 可 能 受到 环境 胁迫 〈Sun etal., 2022) 。 该 比值 的 变化 会 造成 叶 绿 


素 a 对 不 同 波长 吸收 的 改变 ， 进 而 也 会 对 叶 色 产生 影响 ( 李 丽 羡 等 ，2022) 。 
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* 表 示 差 异 显 着 (P<0.05) 。 
* indicates significant differences (P € 0.05) . 


8 YL 和 NYL 的 叶绿素 和 类 胡萝卜 素 含 量 比较 


Fig. 8 Comparison of chlorophyll content and carotenoid pigment content between YL and NYL 


在 叶绿素 合成 途径 中 ， 谷 氮 酸 -让 NA 还 原 酶 (glutamyl-tRNA reductase; HemA) . EK 
啉 原 脱羧 酶 Curoporphyrinogen decarboxylase, HemE? . 、 原 路 啉 原 氧化 酶 (protoporphyrinogen 
oxidase, HemF) 、 原 叶绿素 氧化 还 原 酶 (protochlorophyllide oxidoreductase, POR) 等 9 种 
酶 的 相关 DEGs 都 出 现下 调 表达 , 其 中 编码 HemE 的 基因 下 降 4/5, 编码 POR 的 2 个 基因 表 
达 量 下 调 超过 2/3， 编 码 HemA、HemB C(ó-aminolevulinic acid dehydratase) ~ DVR (divinyl 
protochlorophyllide a 8-vinyl-reductase) 和 ChlI (Mg-chelatase subunit ChlI-1 ) 的 基因 表达 下 


调 都 接近 50% (E 9AD 。 


NS 
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CHN pi ANS 
FXH &W- DFR 
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2 -| 0 1 2 DHQ * Leucocyanidin ————————  Cyanidin 
E à 人 
rysu | ANS 
| | | DHM — E baorik iid » Delphinidin 
YL2 YLA YL3 NYLINYL2 NYL3 DER WW Phinidin /Ns mag 
A. 叶绿素 生物 合成 '， CAO. 叶绿素 酸 酯 a MAN: ChG. 叶绿素 合 酶 ，ChlM. ENIX H hE Fe 
ChI. 镁 离子 整合 酶 亚 基 Chl I-1; CLH. 叶绿素 酶 ; CRDI. 镁 原 中 啉 IX 甲 酯 环 化 酶 ; DVR. 二 乙烯 还 原 酶 ; 
HemB. ;- 氮 基 乙 酰 丙 酸 脱水 酶 ， HemL. 谷 氮 酸 -1- 半 醛 -2,1- 氨 基 变 位 酶 ， B. 类 胡萝卜 素 生物 合成 ; CrtISO. 
类 胡 葛 小 素 异 构 酶 ，LCYB. 番茄 红 素 B- 环 化 酶 ; LCYE. 番茄 红 素 s 环 化 酶 ， PDS. 八 氢 番茄 红 素 脱氧 酶 ; 
PSY. 八 氨 番 茄 红 素 合 酶 ，C. 茶 丙 烷 生物 合成 ， PAL. 茶 丙 氨 酸 解 氨 酶 ，4CL. 4- 香 豆 酸 辅酶 A 连接 酶 ; 
CYP73A. 反 - 肉 桂 酸 4- 单 加 氧 酶 ,， D. 类 黄酮 生物 合成 , ANS. ERKAM: CH. 查 尔 酮 异 构 酶 ; CHS. 查 
尔 酮 合 酶 ; DFR. 二 氧 黄酮 醇 -4- 还 原 酶 ; DHK. 二 氧 山 奈 酚 ; DHM. 二 氧 杨梅 素 ; DHQ. ZARR: F3'H. 
黄 烷 酮 3- 单 加 氧 酶 ，F3'5'"H. 黄 烷 酮 F3'5- 羟 化 酶 ， FLS. 黄酮 醇 合 酶 。 


A. Chlorophyll biosynthesis ; CAO. Chlorophyllide a oxy 


genase ; ChlG. Chlorophyll synthase ; ChIM. 


Mg-protoporphyrin IX O-methyltransferase; ChlI. Mg-chelatase subunit ChlI-1; CLH. chlorophyllase; CRDI. 


Mg-protoporphyrin IX monomethylester cyclase; DVR. divinyl protochlorophyllide a 8-vinyl-reductase; HemB. 


6-Aminolevulinic acid dehydratase ; 


HemL. Glutamate-1-semialdehyde 2,1-aminomutase ; 


B. Carotenoid 


biosynthesis; CrtISO. Carotenoid isomerase; LCYB. lycopene f-cyclase; LCYE. lycopene £-cyclase; PDS. 


Phytoene desaturase ; PSY. Phytoene synthase ; 


C. Phenylpropanoid biosynthesis ; 


PAL. Phenylalanine 


ammonia-lyase; 4CL. 4-Coumarate-CoA ligase; CYP73A. Trans-cinnamate 4-monooxygenase; D. Flavonoid 


biosynthesis; ANS. Anthocyanidin synthase; CHI. Chalcone isomerase ; 


CHS. Chalcone synthase; DFR. 


Dihydroflavonol4-reductase; DHK. Dihydrokaempferol; DHM. Dihydromyricetin; DHQ. Dihydroquercetin; 
F3'H. Flavonoid 3'-monooxygenase; F3'5S'H. Flavonoid 3',5'-hydroxylase; FLS. Flavonol synthase. 
图 9 YL 组 和 NYL 组 叶绿素 、 类 胡萝卜 素 和 类 黄酮 生物 合成 通路 分 析 
Fig.9 Analysis of chlorophyll, carotenoid and flavonoid biosynthesis pathway of in YL and 
NYL 
类 胡萝卜 素 也 是 参与 植物 光合 作用 的 主要 色素 之 一 。YL 中 类 胡萝卜 素 含 量 (0.056 1— 

0.008 9 mgg!) 较 NYL (0.441 5 土 0.083 0 mgg!) 显著 减少 。 在 类 胡萝卜 素 生物 合成 途径 
中 (图 9B) ， 八 氧 番 茄 红 素 合 酶 (phytoene synthase, PSY) 是 该 途径 的 第 1 个 合成 酶 ， 编 
码 该 酶 基因 的 表达 量 下 降 了 4/5， 其 显著 下 调 会 对 后 续 通 路 产生 重要 影响 。 番 茄 红 素 s 环 化 酶 

(Lycopene £-cyclase, LCYE) 是 叶 黄 素 代谢 途径 中 的 关键 酶 ， 众 化 线性 的 番茄 红 素 环 化 形 
成 胡萝卜 素 ， 再 进一步 合成 叶 黄 素 。LCYE 基因 的 表达 量 也 出 现 显著 下 降 ， 下 调 约 4/5， 会 
对 叶 黄 素 生 物 合成 产生 影响 ， 进 而 影响 突变 体 的 叶 色 表现 。 
类 黄酮 合成 的 起 始 代谢 物 为 对 香 豆 酰 辅酶 A 〈p-coumaroyl-CoA) ， 通 过 茶 丙 烷 生 物 合 
成 途径 产生 。 在 YL 组 中 ， 编 码 4- 香 豆 酸 辅酶 A 连接 酶 CA-coumarate-CoA ligase) 的 基因 
4CL 表达 上 调 (图 9C) ， 但 并 未 对 香 豆 栈 辅酶 A 的 含量 形成 显著 影响 。 在 类 黄酮 生物 合成 
初期 阶段 ， 编 码 查 尔 酮 合 酶 (chalcone synthase, CHS) 和 查 尔 酮 异 构 酶 Cchalcone isomerase, 
CHI) 的 基因 表达 显著 上 调 ， 导 致 柚 皮 素 查 尔 酮 Cnaringenin chalcone) 和 柚 皮 素 Cnaringenin ) 
在 YL 中 积累 增多 (FC = 1.05) 。 在 此 后 的 步骤 中 , 许多 关键 酶 的 基因 表达 水 平 也 显著 上 调 ， 
其 中 黄酮 醇 合 酶 (flavonol synthase, FLS) 的 基因 表达 量 上 调 了 6.98 倍 ， 黄 烷 酮 3- 单 加 氧 
H$ (flavonoid 3-monooxygenase, F3'HO 的 基因 表达 上 调 5.77 倍 ， 二 氢 黄 酮 醇 -4- 还 原 酶 
dihydroflavonol4-reductase, DFR) 的 基因 表达 水 平 也 显著 上 调 4.44 倍 。 这 些 基因 的 表达 
上 调 ， 与 YL 中 榭 皮 素 〈quercetin) 、 白 矢 车 菊 素 〈leucoyanidin) 、 杨 梅 素 (myricetin ) 在 
内 的 多 种 黄酮 类 化 合 物 的 积累 有 着 强 相 关 性 (x|>0.8) 〈 图 9D， 图 10) 。 
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CHS. 查 尔 酮 合 酶 编码 基因 ; CHI. 查 尔 酮 异 构 酶 编码 基因 ; DFR. 二 所 黄酮 醇 -4- 还 原 酶 编码 基因 ; FLS. 黄 
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gene; NODS. Naringenin, 2-oxoglutarate 3-dioxygenase encoding gene; LDOX. Leucoanthocyanidin dioxygenase 
encoding gene; UGT. UDP-glycosyltransferase encoding gene. 
10 YL 和 NYL 中 SCMs 与 DEGs 的 互 作 网 络 图 
Fig. 10 Strong correlations between SCMs and DEGs in YL and NYL 


2.32 光合 作用 相关 DEGs 和 SCMs 的 关联 分 析 


在 光合 作用 途径 中 ， 共 有 28 个 编码 PS I. PS ILI LHC 核心 蛋白 的 DEGs 表达 水 平 发 
生变 化 (图 11) 。 其 中 ， 与 PSI 反 应 中 心 相关 的 8 个 DEGs， 都 在 YL 中 表现 下 调 。 与 PS 
开 反应 中 心 相 关 的 10 个 DEGs, 只 有 编码 PS TI sz FoU DI 蛋白 (photosystem II P680 reaction 
center D1 protein ) 的 psbA 基因 上 调 ,$ 个 PS I 放 氧 增强 蛋白 (photosystem II oxygen-evolving 


pur 


enhancer protein) 和 4 个 不 同 分 子 量 蛋白 的 DEGs 都 表现 为 下 调 。 此 外 ， 编 码 叶绿素 a/b 结 
合 蛋白 的 10 个 基因 CLHCA 1, LHCA 2-1. LHCA 2-2. LHCA 4. LHCB 1-1. LHCB 1-2. LHCB 
2. LHCB 4. LHCB 6. LHCB 7) 在 YL 叶片 中 下 调 表 达 。 以 上 结果 表明 ， 下 调 的 与 光合 作 
用 有 关 的 编码 基因 对 YL 叶片 中 叶绿体 的 发 育 有 影响 ， 这 与 YL 组 叶片 叶绿素 含量 降低 的 情 


Ud 


况 相 符 。 叶 绿 素 含量 减少 ， 进 而 引起 YL 组 光合 速率 下 降 〈 表 5) 。 
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图 11 YL 组 和 NYL 组 叶绿素 、 类 胡萝卜 素 和 类 黄酮 生物 合成 通路 分 析 
Fig. 11 Analysis of chlorophyll, carotenoid and flavonoid biosynthesis pathway of in YL and 


NYL 


表 5 YL 和 NYL 的 光合 作用 参数 
Table S$ Photosynthetic parameters of YL and NYL 


EB 


实验 材料 净 光 合 速率 气孔 导 度 胞 间 CO» 浓度 蒸腾 速率 
Materials Net photosynthetic — Stomatal conductance Intercellular CO2 Transpiration rate 
rate (umol-m?- S!) (mmol:m?- S!) concentration (mol:m?- S!) 
(hmolmolD 
NYL 14.475 0.345 324.46 : 3.45 238.67 15.32 7.320125 
YL 6.350 0.210 278.51 +4.62 322.46 +2.18 10.46 +1.42 
3 讨论 


3.1 叶绿素 合成 受阻 和 叶绿体 发 育 异 常 导致 YL 叶 色 黄 化 


叶 色 是 植物 的 一 个 重要 性 状 ， 其 形成 涉及 许多 复杂 的 生物 学 过 程 CWang et al., 2022) 。 
企 高 等 植物 中 ， 叶 色 表 现 主要 决定 于 叶绿素 代谢 过 程 〈 吴 砚 农 等 ，2021) ， 通 常 是 因 基因 突 
变 引 起 的 叶绿素 合成 受阻 或 降解 加 速 ， 导 致 植株 叶 色 黄 化 〈 朱 美玉 等 ，2020) 。 在 YL 突变 


体 的 叶绿素 合成 过 程 中 ， 共 有 HemA、HemB、HemE、ChlI[、POR 等 9 种 酶 的 基因 下 调 表 


达 ， 涉 及 从 工 - 谷 氨 酰 -了 RNA 到 叶绿素 酸 酯 合成 的 大 多 数 过 程 。 其 中 ，HemA EM L-A 
-tRNA CL-glutamatectRNA ) 还 原 ， 是 调控 叶绿素 合成 的 第 一 个 限 速 反应 Zhao et al., 2014)。 


HemB 则 参与 8- 氨基 乙酰 丙 酸 Có-aminolevulinate) 向 胆 色 素 原 转化 ， 也 是 叶绿素 合成 的 关键 


步骤 (Yang et al., 2015) > HemA 和 HemB 基因 在 YL 中 均 下 调 表达 ， 下 调 约 12， 推 测 正 是 


AE 


由 于 该 两 个 前 端 基因 的 显著 下 调 导 致 叶绿素 合成 的 前 体 物质 减少 , 引起 后 续 反应 受阻 , 进而 


使 叶片 整体 变色 。 知 突变 发 生 于 叶绿素 合成 的 后 期 基因 ， 则 突变 体 通常 形成 条 纹 或 斑点 状 


(Sakuraba et al., 2015) 。 例 如 ，HemE 基因 突变 会 导致 玉米 (Zea mays) 叶片 出 现 病 斑 (Hu 


et al., 1998) ，HemF 基因 突变 则 会 引起 玉米 叶片 黄 化 坏死 (Williams et al., 2006) 。 此 外 ， 
CAU 基因 表达 下 调 也 可 能 引起 YL HERM., CH 是 镁 离子 歼 合 酶 (Mgch) 的 三 个 亚 基 之 
一 ，Mgch 则 是 叶绿素 合成 途径 〈 镁 分 支 ) 中 的 第 一 种 酶 ， 也 是 一 种 关键 的 限 速 酶 ( 罗 莎 等 ， 


20150 。 已 有 研究 表明 。Megch 亚 基 编码 基因 ChI, ChD 的 表达 会 受 该 酶 底 物 浓度 的 影响 
(Zhang et al., 2006) > Æ YL 黄 化 突变 体 中 ，Ch11 基因 表达 量 下 调 约 50%， 推 测 可 能 是 由 
TH Bi fk Je MK IX (protoporphyrin IX) 含量 降低 所 致 。 并 且 ， 该 过 程 也 可 能 进一步 抑制 
后 续 反 应 ， 最 终 导致 YL 组 叶绿素 含量 显著 减少 。 此 外 ，YL 中 类 胡 葛 卜 素 含量 的 显著 下 降 ， 


对 突变 体 黄 化 现象 的 形成 也 有 一 定 影响 。 但 是 ， 突 变 体 中 叶绿素 与 类 胡萝卜 素 的 比值 


(Chl/Car) 未 发 生 显著 改变 ， 两 组 比值 分 别 为 8.21 和 8.31， 暗示 Chl/Car 比值 对 YL 的 黄 叶 


表 型 影响 不 大 。 


叶绿体 结构 异常 也 可 能 导致 植株 叶 色 黄 化 〈 江 新 凤 ，2021) 。 在 YL 突变 体 中 ， 有 33 


条 与 光合 作用 相关 的 DEGs 表达 水 平 发 生变 化 ， 其 中 28 条 为 PS I. PS ILI] RCC 和 LHC 的 
编码 基因 。Psa 和 Psb 家 族 中 许多 基因 的 下 调 表 达 , 可 能 导致 PS I 和 PS [中 相关 蛋白 的 功 
能 受阻 。 已 有 研究 证 实 ，PS 了 [和 蛋白 复合 物 的 显著 下 调 会 引起 不 良 基 粒 的 堆 疼 (能 兴 伟 等 ， 
2023) 。LHC 与 色素 结合 起 到 吸收 和 传递 光 能 的 作用 ， 并 能 实现 光 保 护 ， 其 表达 下 降 会 导 
致 叶绿体 中 基 粒 堆积 异常 ， 引 起 黄 叶 表 型 (Kim etal., 2009) 。 在 对 水 稳 黄 叶 突变 体 的 研究 
rH, Wu 等 (2007) 发 现 叶 绿 体 发 育 相 关 基 因 以 及 PS 了 [蛋白 复合 体 基因 的 表达 都 会 受 叶 绿 


素 合 成 中 间 产 物 的 调控 。 在 拟 南齐 研究 中 ，LHCB 基因 的 表达 受到 叶绿素 合成 相关 酶 (如 
Mgch) 的 反馈 调节 (Mochizuki et al., 2001) 。 由 此 推测 ， 由 于 叶绿素 合成 途径 受阻 ， 其 


间 产 物 和 相关 酶 的 含量 发 生变 化 , 从 而 抑制 了 叶绿体 光 系 统 中 RCC 和 LHC 的 基因 表达 , 影 


响 了 类 圳 体 膜 结构 的 形成 ， 导 致 叶绿体 发 育 异 常 ， 以 致 叶 色 黄 化 。 此 外 


u 


， 叶 绿 素 缺 乏 与 叶 绿 


体 发 育 异常 ， 也 导致 YL 突变 体 捕 光 能 力 下 降 ， 对 光合 作用 产生 不 良 影响 。 


3.2 类 黄酮 物质 合成 是 YL 叶 色 黄 化 的 物质 基础 


采用 非 靶 向 代谢 组 (LC-MS) 分 析 ， 在 YL 组 中 共 检 测 到 614 个 SCMs， 主 要 涉及 黄酮 


类 、 氮 基 酸 、 氨 基 糖 或 核 音 酸 糖 等 。 黄 酮 类 物质 中 的 SCMs 数量 最 多 ， 


有 23 个 ， 包 括 白 矢 


车 菊 素 、 杨 梅 素 、 榭 皮 素 等 类 黄酮 化 合 物 及 其 糖苷 衍生 物 。 其 中 ， 吡 喃 酮 啡 肽 A 在 YL R 
变 体 中 上 调 倍数 最 大 ,， 达 2.28 倍 。 该 物质 最 早 发 现 于 黑 加 仑 种 子 提取 物 中 (Lu et al., 2000), 


目前 关于 其 合成 途径 和 形成 机 制 尚 不 清楚 。 黄 酮 类 等 含 氮 化 合 物 中 形成 SCMs 的 积累 , 表明 
YL 突变 体 中 与 碳 和 氮 相 关 的 通路 发 生 了 代谢 重 编 事件。 在 YL 中 ， 参 与 叶绿素 合成 、 碳 固 


定 以 及 光合 作用 的 DEGs 主要 表现 为 下 调 表 达 (HemA. TIM, GAPA, LHC 等 ) ， 会 造成 光 
合作 用 受阻 CX 50 ， 进 而 导致 下 游 的 糖 酵 解 、TCA 循环 等 碳 代谢 过 程 受到 抑制 〈 宋 建 民 等 ， 


1998) 。 糖 醉 解 和 TCA 循环 的 中 间 体 是 黄酮 类 物质 碳 骨 架 的 主要 来 源 
代谢 速率 下 降 会 影响 类 黄酮 的 生物 合成 。 然而， 对 代谢 组 数据 的 分 析 表 


( 刘 健 伟 ，2016) ， 
明 ， 两 个 代谢 通路 的 


中 间 产 物 及 相关 酶 的 含量 在 YL 中 无 显著 变化 。 甚 至 ，YL 组 中 糖 酵 解 途径 的 丙酮 酸 激酶 


(pyruvate kinase) 编码 基因 表达 上 调 了 3.14 倍 ，TCA 循环 中 苹果 酸 合 


酶 Cmalate synthase? 


编码 基因 表达 上 调 达 20.42 倍 ， 瞳 示 黄 化 突变 并 未 引起 叶片 中 糖 酵 解 和 TCA 循环 的 代谢 减 


弱 , 反而 是 受到 促进 。 相同 现象 也 在 赤 霞 珠 葡萄 的 黄 化 突变 体 中 被 发 现 


， 原 因 可 能 是 黄 化 叶 


片 对 碳 源 和 和 氨 源 的 征调 能 力 更 强 ( 陈 迎春 , 2011)， 使 得 黄酮 类 物质 合成 的 碳 骨 架 供应 充足 。 


此 外 ， 由 于 叶绿素 合成 受阻 , 减少 了 对 氨 的 消耗 , 氮 的 积累 也 可 能 在 一 定 程度 上 促进 了 黄酮 


类 物质 的 合成 与 积累 ， 并 为 YL 叶 色 黄 化 表现 的 形成 提供 了 物质 基础 


3.3 转录 因子 参与 YL 叶 色 黄 化 过 程 


江 新 凤 ，2021) 。 


转录 因子 对 生物 的 生长 发 育 具 有 重要 的 调节 作用 ,许多 种 类 还 参与 对 非 生物 胁迫 的 应 答 
(Sun et al., 2022) 。 在 YL 中 ， 转 录 因 子 富 集 最 多 的 DEGs Y MYB, AP2IERF 和 bHLH. 
己 有 研究 表明 ，MYB 家 族 转 录 因 子 能 够 介 导 类 黄酮 合成 途径 中 许多 关键 酶 的 转录 ， 促 进 黄 
酮 类 物质 产生 (An et al., 2017) o bHLH 转录 因子 则 被 证 实 参与 对 环境 胁迫 的 应 答 ， 并 能 协 
同 MYB 调控 类 黄酮 生物 合成 (Liu et al., 2018; Wang et al., 2018) 。 在 了 YL 组 中 , CHS, CHI, 
F3H、FLS 等 类 黄酮 合成 关键 酶 的 基因 表达 上 调 ， 可 能 是 受 MYB 和 bHLH 转录 因子 的 调控 。 
此 外 ，AP2/ERF 家 族 4 个 DEGs 表达 发 生 显著 改 变 ， 该 基因 家 族 主要 参与 非 生物 胁迫 的 应 
激 反 应 ， 如 强 光 照 、 高 温 、 强 光 等 (Wu etal., 2015) . WRKY, HSF 等 与 逆境 胁迫 相关 的 
DEGs 也 出 现 富 集 ， 表 明 YL 黄 化 突变 体 可 能 受到 环境 胁迫 。 进 一 步 分 析 胁 人 迫 和 条件， 可 能 源 
于 强 光 或 水 分 缺失 。 由 于 黄 化 突变 体 缺少 叶绿素 (图 8), 更 容易 遭受 强 光 胁迫 。 江 新 凤 (2021) 
研究 证 实 强 光 胁迫 能 诱导 MYB 基因 上 调 并 参与 植株 黄 化 过 程 。 此 外 ， 更 快 的 蒸腾 速率 也 可 
能 使 YL 遭受 水 分 胁迫 。YL 的 蒸腾 速率 是 NYL 的 1.42 倍 〈 表 $) ， 蒸 膳 速率 过 快 易 导 致 植 
株 缺 水 〈 张 丽 霞 等 ，2021) 。 因 此 ， 推 测 YL 受到 强 光 或 缺 水 胁迫 ， 诱 导 MYB. bHLH 等 转 
录 因 子 表达 上 调 并 参与 类 黄酮 的 生物 合成 ， 促 进 黄酮 类 化 合 物 的 产生 。 


4 结论 


本 研究 通过 代谢 组 和 转录 组 联 用 技术 探索 赤 皮 青冈 叶 色 黄 化 的 形成 机 制 。 结 果 表 明 ， 
YL 金黄 叶 色 可 能 是 受到 叶绿素 合成 受阻 、 叶 绿 体 发 育 异 常 以 及 黄酮 类 物质 合成 加 强 等 因素 
的 综合 作用 。 此 外 ， 在 YL 组 还 发 现 MYB 和 bHLH 家 族 的 基因 表达 水 平 显 著 上 调 ， 证 实 了 
该 两 类 转录 因子 参与 调控 类 黄酮 的 生物 合成 。 本 研究 拓展 了 对 赤 皮 青冈 黄 叶 表 型 形成 机 制 的 
认识 ， 为 更 多 “ 金 叶 ”类 型 园林 植物 资源 的 选 育 工作 提供 理论 依据 。 
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